Enhanced heavy quark-pair production in strong SU(2) color field 
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Non-perturbative charm and bottom quark-pair production is investigated in the early stage ol 
heavy ion collisions. Following our earlier works, the time-dependent study is based on a kinetic 
description of fermion-pair production in strong non-Abelian fields. We introduce time-dependent 
chromo-electric external field with a pulse-like time evolution, which simulates the overlap of two 
colliding heavy ions. The calculations is performed in a SU(2) color model with finite current quark 
masses. Yields of heavy quark-pairs are compared to the ones of light and strange quark-pairs. We 
show that the small inverse duration time of the field pulse determines the efficiency of the quark- 
pair production. Thus we do not see the expected suppression for heavy quark production, as follows 
from the Schwinger formula for constant field, but rather an enhanced heavy quark production at 
ultrarelativistic energies. We convert pulse duration time-dependent results into coUisional energy 
dependence and introduce fiavour-dependent energy string tensions, which can be used in string 
based model calculations at RHIC and LHC energies. 

PACS numbers: 24.85.-|-p,25.75.-q, 12.38.Mh 



I. INTRODUCTION 

The main aim of ultrarelativistic heavy ion collisions 
is to create extreme high energy densities and study the 
deconfinement phase transition of colored quarks and glu- 
ons. Experiments at the BNL Relativistic Heavy Ion Col- 
lider (RHIC) have been investigated the center of mass 
colliding region up to ^/s = 200 AGeV and detectors 
at CERN Large Hadronic Collider (LHC) are ready to 
explore the energy range to y/s = 5500 AGeV. At such 
high energies the colliding nuclei are two colliding sheets 
of nucleons with a huge Lorentz-contraction (7cm = 100 
at RHIC and = 2750 at LHC), surrounded by a 
gluon cloud. Their overlap results in a strong chromo- 
electric and chromo-magnetic field to be built up. Par- 
ticles, namely gluons and quark-antiquark pairs are pro- 
duced from this strong field, similarly to the Schwinger 
mechanism in quantum electrodynamics (QED) The 
particle production rate depends on the field strength, 
which is varying in time. The soft particles produced 
in such a non-perturbative way form the bulk of the 
wanted quark-gluon plasma, after they successfully ther- 
malized. Light and strange quarks loose most of their 
original properties during thermalization, but charm and 
bottom quarks conserve certain characteristic properties, 
which can be studied after the whole evolution of the 
heavy ion collisions. In this paper we investigate the pri- 
mordial non-perturbative production of heavy charm and 
bottom quarks, and explore the early stage of heavy ion 
collisions. 



Recently the study of heavy quark production has re- 
ceived wide interest, because open charm has been mea- 
sured at RHIC in d-hAuJlL Cu-hCu [1], and Au-hAu [i| 
collisions. The STAR [ill, H and PHENIX [1, @ ex- 
periments have obtained different results (within a fac- 
tor of 3), which finding opened vivid experimental and 
theoretical discussions [tJ. A review on heavy-fiavour 
production has been pubfished recently Theoreti- 
cal calculations based on perturbative quantum chromo- 
dynamics (pQCD) at fixed order next to-leading logar- 
itms (FONLL) have found that the measured total charm 
cross section only comparable with the upper limit of 
the FONLL calculations. Thus there is a room for non- 
perturbative production channels. However, the situa- 
tion is very complex, as it was discussed in a recent pa- 
per 0: measured heavy quark radiative and coUisional 
energy loss in heavy ion collisions must fit into the pic- 
ture, as well as the regeneration of charm hadrons in 
quark coalescence channels. The interest on this question 
supports the importance of our investigation of primor- 
dial non-perturbative heavy quark production. 

Theoretical descriptions of particle production in high 
energy pp collisions are based on the introduction of chro- 
moelectric fiux tube ('string') models 0, [O, d, [H, [ll ■ 
String picture is a good example of how to convert the ki- 
netic energy of a collision into field energy, than later on 
gain the stored kinetic energy back. However, at RHIC 
and LHC energies the string density is expected to be so 
large that a strong collective gluon field will be formed 
in the whole available transverse volume. Furthermore, 
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the gluon number will be so high that a classical gluon 
field as the expectation value of the quantum field can 
be considered in the reaction volume [isl . 

Al- 

ternatively at extremely high energies, nucleus nucleus 
collisions can be described as two colliding sheets of Col- 
ored Glass Condensate. In the framework of this model it 
was shown that in the early stage of collision longitudinal 
color-electric and color-magnetic fields are created p7| . 
The properties of such non-Abelian classical fields and 
details of gluon production were studied very intensively 
during the last years, especially asymptotic solutions (for 
summaries see Refs . fisl. I20I). Fermion production was 
calculated recently [llnillii]. 

Fermion pair production together with boson pair pro- 
duction were investigated by different models of particle 
production from strong Abelian 24i.i25i. j26.i27u i28i.i2S.i3Q. 
|3ll, m m, ill and non-Abelian [sl, S [33| fields. These 
calculations concentrated mostly on the bulk properties 
of the gluon and quark matter, the time evolution of the 
system, the time dependence of energy and particle num- 
ber densities, and the appearance of fast thermalization. 

In our previous papers (see Ref. [H, [s^) we investi- 
gated massless fermion and boson production in strong 
Abelian and non-Abelian external electric field. During 
these calculations we have realized, that the role of mass 
becomes unimportant when the coUisional energy is in- 
creasing and the the pulse duration time becomes com- 
parable to the inverse quark mass ,40]. In this paper 
we describe strange, charm, and bottom quark-pair pro- 
duction. Motivated by the problems raised in Ref. [9] 
we investigate the role of string tension in the Schwinger 
mechanism for heavy quark pair production. 

The energy dependece of the string tension was inves- 
tigated earlier [l^. Instead of the usual k ~ 1 GeV/fm 
value the much higher effective string tension, k ~ 5 — 12 
GeV/fm appeared in the calculations. This open ques- 
tion motivates our investigation also. 

In this paper we solve the kinetic model in the pres- 
ence of an external SU(2) non-Abelian color field. We 
focus on particle production with finite mass at differ- 
ent duration time of the quickly changing external field. 
Section 2 summarizes theoretical background of the ki- 
netic equation for color Wigner function. In Section 3 
we consider the kinetic equation for pure longitudinal 
time-dependent SU(2) color field, which will be solved 
numerically. In Section 4 we summarize and discuss our 
numerical results from general point of view. We discuss 
the collision energy and pulse duration time-dependence 
of heavy quark production. In a phenomenological way 
we introduce flavour specific string tensions to connect 
our numerical results and Schwinger estimate. 

In Appendix A we derive general kinetic equation for 
the Wigner function starting from QCD Lagrangian. In 
Appendix B we show exact solution of the kinetic equa- 
tion for SU(2) color case, that supports our numerical 
calculation. 



II. THE KINETIC EQUATION FOR THE 
WIGNER FUNCTION 



The equation of motion for color Wigner function 
VF(k, t) in the gradient approximation reads (see Ap- 
pendix A for details) 

= zA:,{7°y , W} - ^m[7^ W] + ig [A, , [7°7\ W]] . (1) 

Here m denotes the current mass of the fermions, g is the 
coupling constant, is the 4-potential of an external 
space-homogeneous color field and Fi^„ is corresponding 
field tensor 



F^u = d^A^ - d^A^ - ig[A^, A^]. 



(2) 



The validity of gradient approximation requires that the 
Wigner function is sufficiently smooth in momentum 
space and the field strength varies slowly in coordinate 
space. The corresponding characteristic lengths must 
satisfy the following relation {Ap)w{^x)f ^ 1, where 
{Ap)]Y is connected to the momentum gradients of the 
Wigner function and {Ax)f to space gradients of the 
field. 

The color decomposition of the Wigner function with 
SU(A'c) generators in the fundamental representation is 
given by 



1,2, 



1 



(3) 



where is the color singlet part and W"" is the color 
multiplet components. It is also convenient to perform 
spinor decomposition separating scalar a, vector 6^, ten- 
sor Cpjy, axial vector and pseudo-scalar parts e: 



+ bl}"!^ + c'Jy^" + dl^''-f>'-f^ + ie"l"7^ (4) 



The asymmetric tensor components of the Wigner func- 
tion is convenient to decompose into axial and polar vec- 
tors cj = c'^ and C2 = ^e°"''"'C(^p correspondingly. 



III. KINETIC EQUATION IN SU(2) WITH 
COLOR ISOTROPIC EXTERNAL FIELD 

After the color and spinor decomposition of the equa- 
tions for the Wigner function in case of pure longitudi- 
nal external SU(2) color field with fixed color direction 
A"^ = Aln", where n"n" = 3 and dtn'' = [111, we 
obtain the following system of equations for the singlet 
component 



* ^ 4 ^dk. 



-4kc2 - 2m(ig, 



(5) 
(6) 
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dth' + '^Et^h'> = 2[k X d^] - 
4 dkz 

dtdl + ^-K^dl = 2me^ 



35 po d 



and the triplet components 



d 



dte^ + gEt^^e^ = -4k4 - 2mdg, (14) 

5ib* + gEt^h' = 2[k X d*] + 4mc?, (16) 

9td^ + g£;:^dg^''^-2me^ 

a*do + .9£;:Ad^^2[kxb*], 

dtc\ + gEt^c' = a*k - mb^ 
dfcz 



(7) 




-^El—b\ = ~2Kdl^Amc\, 


(26) 


(8) 


dtdl^ 


- = 2(fc.6i - k^bw 


(27) 


(9) 


dtcl + 




(28) 


(10) 






(29) 


(11) 


and for the triplet components 




(12) 




fg£;:^a^ = -4(^,4 + fc^cl), 


(30) 




dthl + gEt-^hl = 2k^dl + Amcl, 


(31) 


(13) 


dtbl 


+ 5i?:^fol = -2A;,rf^+4mcl, 


(32) 


(14) 


dtdl 


^gEt^dl^2{k,b\-k^bl), 


(33) 


(15) 




i-gEt-^^cl^aOk^-mbt, 


(34) 


(16) 




+ gEt-^cl^a^k^-mbl ; 


(35) 



dt4+gEt-^^c^,=e-k 



(17) 
(18) 
(19) 
(20) 



The distribution function for massive fermions is com- 
pletely defined by components a,b [s^: 

, n. .^ _ ma-(k,t)+kb-(k,t) , 1 

~ ^k^ + 2 ' ^^^^ 



where ^(k) — Vk^ + m^. Thus for time- and 
momentum-dependent distribution functions scalar a, 
vector 6^, axial vector d^, and tensor c components of 
the Wigner function are needed, only. 

The initial conditions for the Wigner function in vac- 
uum reads (see Appendix A) 



1 m 

" = -2^^ 

b^ ^ l^. 

2uj 



(22) 
(23) 



Considering vacuum initial condition symmetry we ob- 
tain the following equations for the singlet (we redefined 
c = Ci to simplify reading) 

dta' + ^Et^a'' = -A{k,4 + k^cl), (24) 
4 okz 

dtbl + ^El-^bl = 2k^d% + 4mcJ, (25) 



where we introduced the following vector decomposition 



n X 



(36) 



Here the unit vector coUinear to the field direction, n, is 
given by n = E^/|E*| = (0,0,1) and k^ = (fci,/c2,0). 
As follows from Eqs. (|24ll35p . the axial part of vector bx: 
tensor components Cx , longitudinal d^ and perpendicular 
d^ parts of axial vector components do not contribute to 
the evolution of the distribution function. 

These equations could be further simplified and trans- 
formed to ones that are similar to Abelian case (see Ap- 
pendix B for details). 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. General results 

In Ref. [11] we have solved the above equations for 
massless (light) quarks and described their longitudinal 
and transverse momentum distributions. Here we focus 
on the integrated particle yields and discuss the obtained 
results, focusing on massive (heavy) quark production. 

In the numerical calculation we have used the fol- 
lowing parameters: the maximal magnitude of the field 
Eo = 0.68 GeV/fm; the strong coupling constant g — 2; 
the current quark masses m„^d = 8 MeV, — 150 MeV, 
TOc = 1-2 GeV, mi, = 4.2 GeV for light, strange, charm 
and bottom quarks, correspondingly. The value of max- 
imal magnitude of the field corresponds to the effective 



4 




FIG. 1: The dimensionless quark number density for differ- 
ent flavours, nq{t)Eg as a function of dimensionless time 



tE^^ for different pulse duration, tEq'^ = 0.1 (solid lines) 



1/2 



and tE, 



1/2 



: 0.5 (dashed lines). 



string tension k ~ 1.17 GeV/fm. The reason why we use 
this value will be explained below, in Subsection IIVBI 

The particle production is ignited bya pulse-like color 
field simulating a heavy ion collision [39| : 



E^it) ^Eo- [l-tanh2(Vr)] , 



(37) 



where r is a pulse duration time. 

The suppression factor of heavier quark Q to light one 
u is defined in the asymptotic future (c.f. [1]), t ^ r, as 



7^ = lim nQ{t)/nu{t). 



(38) 



Here Uq (t) is the number density of corresponding quarks 
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FIG. 2: The dimensionless quark number densities at the 
asymptotic final state, nq{t ^ ''")-Bg as a function of pulse 

1/2 

duration time tEq' . 
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FIG. 3: The suppression factor, 7*^, at different pulse dura- 
tion time as a function of the dimensionless mass parameter, 
m/Eo' . 



given by 



(Pk 



(39) 



where q denotes different quark fiavours, q=u, d, s, c, b. 

In Fig. [1] the time evolution of quark number densi- 
ties, Hq, are displayed for different pulse duration time, 
1/2 

tEq = 0.1 and 0.5. For short pulse the quark number 
densities are comparable with each other (solid lines). 
In this case the particle production happens during the 
whole evolution of the field. In contrast to this, for 
long pulse, the number of produced charm and bottom 
quarks becomes negligible in the final state, because their 
production is balanced by annihilation. In U(l) color 
case the annihilation term can be identified clearly, see 
Ref. [13. 

This dependence on the pulse duration time is also 
demonstrated in Fig. [51 This figure clearly displays that 




m X 



FIG. 4: The suppression factor for heavy flavours at different 
values of mr. 
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both charm and bottom quark production is substantially 
enhanced in the cases of short pulse. This enhancement 
has a maximum at r ~ 0.1 ^/Eq for charm and at even 
smaller r for heavier bottom quark. In opposite to the 
heavy quarks, light and strange quark productions are 
increasing with the pulse duration time without any local 
maximum. 

Fig. [2] displays an important result that at small value 
of T the heavy quark production has a pronounced max- 
imum with value well beyond the well known asymptotic 
Schwinger estimate, see Ref. 

We further investigate the suppression factor and its 
dependence on pulse duration time and quark masses. 
Fig. [3] displays the suppression factor, 7*3, for different 
pulse duration time, as a function of quark mass. For 
short pulse the suppression factor is decreasing almost 
linearly with increasing quark mass value. For long pulse 
we can see a very fast exp{— m^/_Bo}) drop, which is 
consistent with the Schwinger formula. 

In Fig. m we demonstrate the pulse duration time de- 
pendence of the suppression factor for different flavours. 
As it can be seen the dramatic change of the suppres- 
sion factors for heavy quarks happens in the region 
0.1 < mr < 1. Actually, as it is shown in Appendix 
B, there are two dimensionless parameters that control 
the behaviour of the particle production. One is the di- 
mensionless pulse duration, mr; the other is adiabaticity 
parameter, F^ ~ m/{EoT). The Schwinger formula is 
valid for the combination mr ^ 1 and Tk ^ 1- 

Note, that at short pulse, the relative charm and bot- 
tom production is surprisingly large, which does not fol- 
low any earlier expectation. 

B. Energy and pulse duration time dependence 

To fix free parameters we use the following simple 
model. The field depends on two unknown parameters, 
T and Eo. We will fix them as the best fit of the sup- 
pression factors for primordial strange and charm quarks 
obtained in a quark coalescence calculation [U, [l^l at 
RHIC energy, ^/s — 200AGeV. The suppression factors 
are 7^ = 0.88 and 7^= = 6 • lO'^. The best fit reads 
Eo = 0.68 GeV/fm and tq = 0.134 fm/c. Surprisingly, 
our simple model provides reasonable values for these pa- 
rameters. 

From intuitive reasons the duration of field pulse is 
proportional to the time of two Lorenz-contracted heavy 
ions pass each other at almost speed of light, i.e. 

2^ 

T~a , (40) 

Tcm 

where R is the radius of a nuclei, 7cm — ^/s/{2GeV) is 
the gamma-factor, a is an unknown proportionality coef- 
ficient. In case of gold-gold collision at RHIC energy we 
obtained a — 0.96 from the best fit values given above. 
We further assume that a weakly depends on the colli- 
sion energy and this dependence can be neglected. Thus 



s [A GeV] 




0%l 0.01 0.1 

T [fm/c] 

FIG. 5: The suppression factor, 7*^, for strange and heavy 
flavours as a function of pulse duration time and collision 
energy. The number marking the curves for strangeness sup- 
pression correspond to different values of /3, see Eq. (|4ip 



having the value of a in hand we can transform the dura- 
tion time dependence to the collision energy dependence. 
Although this conversion is oversimplified (e.g. it does 
not take into account any stopping effects), we expect to 
obtain results of the right order of magnitude. The ex- 
tracted numbers make possible to interpret our numerical 
results on realistic basis, namely energy scales. 

In Fig. [5] we show the time and energy dependence 
of the suppression factor for strange, charm and bottom 
quarks. To demonstrate how robust our results we cal- 
culated the suppression factors assuming different pulse 
duration time (corresponded to different collision ener- 
gies) dependence of the field pulse magnitude Eq (t) . We 
consider three special cases 

Eoir) = Eo ■ (^^y , /3 = 0, 1/2, 1. (41) 

The constant Eo{t) is recovered if /9 — 0. The second 
choice, P = 1/2, corresponds to finite number of quarks 
for short pulses, r ^ 0, (see Appendix B, Eq. ([55])). 
Finally, /3 = 1 results in divergent number of quarks for 
r ^ (see Appendix B, Eq. (PT|) ). 

In Fig. [5] we include these numbers, (3, to distinguish 
between different cases for the strange quark suppression 
factor. For the heavy quarks all three curves are indis- 
tinguishable. As it can be seen, for strange quark the 
difference between the above cases is only important for 
low energy collisions, or long pulse duration. At constant 
Eq, and in the limit r — > 00 the suppression factor for the 
strange quark production tends to the Schwinger limit, 
~ 0.74. For charm and bottom quarks the Schwinger 
limit is negligibly small, see the suppression factor, 7*5, 
for experimentally favoured energies in Table I. 
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7oo 


130A GeV 


200A GeV 


lA TeV 


2A TeV 


5. 5 A TeV 


s 


0.74 


0.84 


0.88 


0.96 


0.98 


0.99 


c 


3-10"^ 


9-10"^ 


0.06 


0.66 


0.82 


0.91 


b 


~ 


~0 


10"" 


0.15 


0.45 


0.72 



TABLE I: The suppression factor, 7^, for experimentally 
favoured energies. The calculations are done with a string 
tension k ~ 1.17, /3 — 0. The result obtained by Schwinger 
formula is denoted by 700. 



C. Effective string tension 

In the previous subsections we have obtained that 
the time dependence of the strong color field can in- 
duce an enhancement of heavy quark production (es- 
pecially in a case of pulse duration times corresponded 
to high energy collisions). In usual string-based mod- 
els [10, [m, m, m, these abundances can be intro- 
duced by increasing the effective string tension in the 
Schwinger formula. However our results show that dif- 
ferent effective string tensions should be introduced for 
different quark flavours. Here we calculate these string 
tension values and offer them for later use. 

Let us recall the Schwinger formula for particles with 
mass m 



dN _ 
dt(Px ~ 47r3 



exp 



(42) 



Here k, is the string tension. According to this formula 
the suppression factor of heavier [Q) to light quarks (q) 
is given by 



exp 



(43) 



This formula is valid for the case of arbitrary N in SU(N), 
see Ref. [H. 
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FIG. 6: The effective string tension, K,^ff, defined by Eq. (|44|) 
as a function of pulse duration time, r, and collision energy. 



At first we extract an effective string tension assuming 
its common value for light and heavy quarks. Provid- 
ing our numerical calculation of suppression factor (see 
Fig. O for Q-fiavour, 7'^(t), we solve the equation 



7c 



(44) 



to find the effective string tension dependence on pulse 
duration time, r, and subsequently on the collision en- 
ergy, -y/s, for given quark flavour. 

Fig. [6] displays results in a logarithmic scale. As it 
can be seen, the values of are very much different for 
strange, charm and bottom quarks. In Fig.[7]we show our 
results in a linear scale focusing on RHIC energy range. 
These large string tension values should be applied to 
light quark production as well, but we see three different 
values. This indicates the need of introduction of effective 
string tension in a different way. 

We keep the usual string tension for light quark ~ 
1.17 GeV/fm and introduce "flavour specific" effective 
string tensions for heavier fiavours. On the basis of 
Eq.(|321) we obtain the suppression factor 
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exp 



(45) 



"off , 



Extracting such a "fiavour specific" effective string ten- 
sions from the numerically calculated values we obtain 
Fig. [5] and Fig. [51 We would like to emphasize the dif- 
ference between n'^^ and k^, which is demonstrated by 
their different values. For easier use we generated the 
Table II displaying the "ffavour specific" effective string 
tensions for experimentally favoured collision energies. 

The values for strange quark are approximately energy 
independent. However for charm and bottom quarks we 
get large values (^ 6 — 38 GcV/fm). Earlier analysis [l^ 
indicated the evidences of large values for string tensions, 
K ~ 5 — 12 GeV/fm. The applicability of such large val- 
ues can be verified after performing proper string model 
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FIG. 7: The same as in Fig. [6] in linear scales, but zoomed 
to RHIC energy range. 
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FIG. 8: The "flavour specific" effective string tension, k'^ff 
as a function of pulse duration time, r, and collision energy. 





130A GeV 


200A GeV 


lA TeV 


2A TeV 


5. 5 A TeV 


u, d 


1.17 


1.17 


1.17 


1.17 


1.17 


s 


1.24 


1.26 


1.32 


1.33 


1.34 


c 


3.32 


4.2 


6.1 


6.3 


6.5 


b 


10.3 


14.7 


32 


36 


38 



TABLE II: The "flavour specific" effective string tension, k^g, 
values (in GeV/fm) for experimentally favoured energies. 



based calculations. As it was demonstrated in Ref. 
the available experimental data on charm production at 
RHIC (for y/s = 200A GeV) is successfully described by 
an effective string tension ~ 3 GeV/fm (the correspond- 
ing value in Table II is close to this result). 



V. CONCLUSION 

We have calculated non-perturbative quark pair pro- 
duction in time-dependent strong non-Abelian SU(2) 
fields. Applying a pulse like time evolution and investi- 
gating the influence of pulse duration time, we observed 
that light and strange quark-pairs are produced as we 
expected, approaching the Schwinger limit. Charm and 
bottom quarks followed this behaviour for long pulses. 
However, for short pulses we did not see the expected 
heavy quark suppression, connected to the large quark 
mass in the Schwinger estimate. Indeed, the large value 
of inverse pulse duration time, overwhelming the mass of 
the heavy quark, l/r 3> rric, determines the quark-pair 
production. We obtained enhanced heavy quark produc- 
tion at small duration time of the pulse, which can be 
connected to ultrarelativistic heavy ion collisions. 

On the basis of our numerical results obtained from 
the kinetic equations we defined "flavour specific" effec- 
tive string tension values to describe the enhanced heavy 
quark production. The validity of the obtained values 
must be verified by string based Monte-Carlo calcula- 
tions. However our values seem to be reasonable in com- 
parison to previously published string model results. 

Finally, we would like to emphasis the strength of our 
model demonstrated in this paper: we are able to de- 
scribe non-perturbative particle production (in strong 
non-Abelian field) in a wide energy range, simulating the 
environment of heavy ion collisions at different energies. 
The flexibility of our model is very much favoured to un- 
derstand the experimental data and the physics behind 
them at different coUisional energies. In the widely used 
CGC model jl^ asymptotic solution could have been ex- 
tracted displaying heavy ion reactions with infinite col- 
lisional energy, only. Our results, the obtained energy- 
dependent enhancement of non-perturbative heavy quark 
production, display the complexity of strong field physics 
and the importance of continuously varying energy de- 
pendence. 
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FIG. 9: The same as in Fig. |8]in linear scales, but zoomed 
to RHIC energy range. 
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Appendix A. Kinetic equation for the Wigner 
function 



The derivation of the kinetic equation for the Wigner 
function in non-Abelian case was discussed in details in 
Ref. [3^ . The covariant proper time formulation of the 
kinetic equation was extended from earlier Abelian ver- 
sion [i^l to non-Abelian one [s^. In our investigation. 
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here and earlier [3^, |40| , we use this non- Abehan version. 
However, the original paper of Prozorkevich et al. [s^ 
contains a few misprints in important equations, which 
may confuse the reader and questioning the validity of 
our work. To avoid this confusion wc shortly display the 
most important steps for derivation of the kinetic equa- 
tion for the Wigner function in non-Abelian case. 

The starting point is the equations of motion for the 
quark field operators, obtained by the variation of the 
QCD Lagrangian. The last reads 



(46) 



where -D^ = — igA^ is the covariant derivative, = 
A^t" is the 4-potential of color field, m is the current 
quark mass and g is the QCD coupling constant. The 
field tensor F^^^ is given by 

F^, = [D^, D,] = d^A, - d,A^ - ig[A^, A,]. (47) 

The generators of SU(N) group in the fundamental repre- 
sentation ta = Aa/2 are expressed in terms of Gell-Mann 
matrices for SU(3) and Pauli matrices for SU(2). 

The wanted equations of motion for spinor field oper- 
ators read 

= --/°YD,'ip-im-/Qijj + igAo'tp, (48) 
dt'ip = 4>j°j'-D* + im^p-fo - ig^'Ao. (49) 

Here the covariant derivative D* in the second equation 
acts to the left. 

The single-time Wigner function ^36| is defined by 



W 



(50) 



J d^^ye^P^U (x, X + I) p(x, y, t)U (x 



the unitary link operator C/(xi, X2) is introduced to main- 
tain gauge invariance of the Wigner function. The link 
operator is given by 



U{xi,X2) = ex.p ^9 J dz^'A^'iz) 
The oneparticle density matrix reads 



ab 
Pik 



(51) 



(52) 



We apply the first derivative w.r.t. time to both sides 
of Eq. ((50|) and take into account that the variation of 
the link operator is given by 



5U{xi,X2) = 
igSx'^A^{xi)U{xi,X2) 



igSx2U{xi,X2)Af^{x2) - 



dsU{xi, z)F^i,{z)U{z,X2){x'^ - x'^) x 



Defining the Schwinger string 

A^-^\X2) = U{XI,X2)A{X2)U{X2,XI) (54) 

we rewrite Eq. (j53p 

5U{xi,x2) = [igSx'^^A^.ixi) - ig5x'iA^^'\x2)- 
-ig f dsF}^j\z){x'i - 4') [6x^2 + <5x'i - fa^)]} X 

U{xi,X2). 

Keeping in mind this variation we obtain the wanted 
equation for the Wigner function: 



dtW + 



+ im[-fo, W] - iki {"fol\ W} + 

\-^j\s{Ef\.-)W + WEf\.+ )+ 

\Ft\^-) ([1^,7V] - s {T^,7V}) (55) 
+ ^([T^,7V]+5{W^,7V})i^i:'(z+) 



-ig[A^,W]--igWAA,,W]]^Q. 

Here we used the equation of motions for fermion field 
operators and the following identity 

j /(y) cxp (iky) d^y ^ J f (~*^) '^^P (^^^^ 

(56) 

that is valid for any analytical function /(y). The argu- 
ment of field tensor is given by 



5* = x±i 



s d 
'29k' 



(57) 



As we can see in Eq. (j55p . the obtained kinetic equation 
is nonlocal and difficult to solve numerically. With help 
of gradient expansion we can derive local approximations. 
In our calculations we use space homogeneous fields. In 
this case the kinetic equation take the form of Eq. ([T]), 
that we used in this paper and previous ones [H, |40] . 

The Wigner equation for the vacuum state follows from 
the definition displayed in Eq. ((50|) . Indeed, using usual 
anticommutation relations for fermion field operators we 
obtain the Wigner function in vacuum 



ab _ ™ + k7^a6 



2uj 



(58) 



{8x2 + s{Sxi — ^2:2)} , z — X2 + s{xi — X2). (53) 



Appendix B. Exact solutions for SU(2)-color case 

In Section HI the kinetic equation for the Wigner func- 
tion was solved numerically. However, after taking into 
account additional symmetries of the external field we 
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can discover further simplifications and even obtain ex- 
act solutions. 



To demonstrate this fact we rewrite Eqs. (|24II35|) ex- 
plicitly for the distribution function /(k, t): 



dtf 

dtv' 

dtu" 

dtf 

dtv^ 



d 



dtu'^ + gEt—u^ 



Here the following new functions were defined 



r = 



W = 



gEei 



kz I ma" 

m 
— I 

e± 



k I b", 



(59) 
(60) 
(61) 
(62) 
(63) 
(64) 

(65) 
(66) 
(67) 
(68) 



The naming scheme is chosen to be consistent with the 
U(l) case of our previous work 38]. The transverse one- 
particle energy in Eq. ((68|) is defined by e±_ ~ •\//c'[ -I- m^. 

As it follows from Eq. ((23|) the vacuum state corre- 
sponds to zero initial conditions for the functions /, 



Note, that Eqs. (|59II64|) has the same number of equa- 
tion for massive and massless particles. Thus the mass- 
less limit does not lead to any further simplifications. 
The system (|59ll64p can be transformed to more conven- 
tional form, allowing the solution on characteristics. For 
that we introduce the following new functions 



U ± U 

2 



The equations for these functions read 



d_ 

dkl 



2^ ^dkz 



-2ujU^, 



2uV^ 



(69) 
(70) 
(71) 

(72) 

(73) 
(74) 



The equations for "(-I-)" and "(-)" functions are com- 
pletely factorized and can be solved independently. The 
distribution function / is obtained as 



/ 



F+ +F- 



(75) 



The equations (|72ll74p are very similar to those obtained 
in the Abelian case (see Eqs. (76-77) in [s^l and Eqs. 
(22-24) in 38]). However, only / is physical quantity, 
while the functions F^ carry intermediate information. 
Nevertheless, from mathematical point of view there is no 
difference and this analogy allows to exploit U(l) solution 
to obtain exact analytical results for SU(2)-color case, as 
we demonstrate below. 

In the case of a time reversal symmetry of the field 
strength, E{t) = E{—t), the corresponding components 
are equal to each other, e.g. F"*" = F~. Thus the distri- 
bution function is defined as / = F^ = F^ . Furthermore 
in the Abelian case it is known that for the field with form 
of Eq. (|37p an analytic solution of Dirac equation exists, 
as well as analytic solution of the corresponding kinetic 
equation (see e.g. Ref. [i^). Following this analogy we 
can obtain an analytic solution of Eqs. (|72j|74p in the 
asymptotic state t ^ t: 



f 



sinh {Tr{9 - /i+ -t- ^J,~)) sinh {tt{0 + ^J,'^ - 
sinh(27r/i+) sinh(27r/i^) 



Here we introduced the following notations 



(76) 



]j(^k,±^m^T^ +kl + m^, (77) 



Eq 2 2 
—niT, 



Er 



2m^ 



(78) 
(79) 



The Schwinger limit can be readily obtained from 
Eq. ([7S|) . Indeed, for pulses longer than any scale in 
the system, r ^ maxjE'p ^^^,m~^/^}, we can use the 
expansions 



± (m+ - 



Eo 



± ^0 2 2 _i_ , 

— m T ± ksT 

2Efr 



^jnV^k,r + oy-), (80) 



%/3 4g£;o 



In the leading order of expansion parameter r we obtain 
the following distribution function 



/ ~ exp -TT 



2ei 



V^gEo 



(82) 



After the integration w.r.t. momentum we obtain SU(2) 
version of the Schwinger formula 



16n 



o)~ I 27rm 
exp 



VigEoJ ' 



(83) 
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where we have taken into account the replacement 
/ dkz V^gEoT/2. Thus in SU(2) the suppression fac- 
tor of heavy particles with mass mg to light particles 
with mass ruq in Schwingcr limit is given by 



exp 



(84) 



The analytical result for the number density of quarks 
can be also obtained if the following inequality is satisfied 



Ecr 1 r T 

Tk = — > maxjl, mrj, 

E mr 



(85) 



where we introduce the adiabaticity parameter [46 
In this limit the dictribution function is given by 



2lj 



-gEr^csch (ttwt) 



(86) 



The momentum integration of this distribution function 
can be done analytically in the following cases: 

a) Long pulse duration and undercritical field. For 
long pulse duration the constraint mr 3> 1 in 
Eq. (I85p is satisfied for undercritical field, E/ Ecr ^ 
1. After expanding the distribution function in 
Eq. (|86p and performing momentum integration, 
we obtain the number density 



n = 



3(mT)3/2 



T{gEo)'^ exp(— 27rmr). 



(87) 



The condition of long pulse duration and undercrit- 
ical field could be realized for the collision of heavy 
ions in SPS energy range (assuming that the phys- 
ical picture of classical gluon field is still valid) for 
light quarks, or for charm and bottom quarks at 
and below the RHIC energies (the string tension is 
the order of 1 GeV/fm.) 

b) Short pulse duration. In the opposite limit, mr ^ 
1, we obtain the number density 



For the parameter set we used in the main part of 
the manuscript the field magnitude, Eq, is about 
five times higher than the critical one for the 
strange quark. The number of produced strange 
quarks follows the Eq. (|88p for pulse duration at 
least five times less then inverse mass of strange 
quark t <C (Sto^)"^ ~ 0.3 fm/c As we estimated 
in Eq. (HDl) the pulse duration time for RHIC is 
about 0.1 fm/c, that is only three times less than 
(5ms)~^. Thus the expression ([55)) is vahd for the 
strange quark only at higher than RHIC energies. 

For heavier particles, e.g. charm quark, the require- 
ment rricT <C 1 leads to smaller values of r. Indeed, 
we can rewrite 



rricT = — X (jTisT) <C 1. 



(89) 



This expression shows that Eq. ((55|) becomes valid 
for charm quark on a shorter scale of pulse duration 

time {rrLc/ms ~ 8). 

For light quarks the condition m„T ^ 1 is trivially 
satisfied at RHIC energies. However now the condi- 
tion (|85p plays more important role. From Eq. (j85|) 
we obtain an estimate for the validity of Eq. 



T < 4 • IQ-^im/c. 



(90) 



In this limiting case the strange suppression factor 
tends to unity, that supports our numerical result 
in Fig. [S]. 

One more analytical result can be obtained from the 
general solution Eq. ([75]) . If the duration time of the 
pulse tends to zero, but amplitude is increasing as Eq = 
Aqt~^ {Aq is constant), then the distribution function is 
given by 



If LO 



3 42 
4^0 



(91) 



127r 



(88) 



In this case the number density of produced parti- 
cles depends linearly on the duration time, r, and 
is independent of the particle mass, m. 



where lo± = 2^±/t. Since the momentum integral from 
the above distribution function is divergent, this approx- 
imation results in infinitely many new quark-antiquark 
pairs in unit volume. This is understandable since we 
pump infinite energy to the system in this special case. 
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